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Abstract 



The decay rate of the neutral long-lived K meson into the e + e~e + e~ final 
state has been measured with the NA48 detector at the CERN SPS. Using 
data collected in 1999, a total of 132 events has been observed with negligible 
background. The total number of kaons was determined to be 5.1 x 10 10 . 
This observation corresponds to a preliminary branching ratio of T(Kl — > 
e + e-e+e-)/r(K L -» all) = (3.67 ± 0.32 stat ± 0.23 sys ± 0.08 norm ) x 10" 8 , 
where the first error is statistical, the second systematic and the third error 
is due to the uncertainty in the normalization. 



I. INTRODUCTION 

The Kl — > e + e~e + e~ decay is expected to proceed mainly via the intermediate state 
Kl — > 7*7* [|T]-§| and thus depends on the structure of the Kl — > 7*7* vertex. Phenomeno- 
logical models include vector meson dominance of the photon propagator QCD inspired 
models ||, intermediate pseudoscalar and vector mesons || and models based on chiral 
perturbation theory j7|. The probability for both virtual photons to convert into e + e~ pairs 
is calculated to be in the range (5.89 - 6.50) x lCT 5 [|H. The chiral model prediction of § 
corresponds to T(Kl — > e + e~e + e~)/T(KL — > all) = 3.85 x 10~ 8 , including the effect of a 
form factor, which increases the width by 4%. The interference term due to the identity of 
particles has been calculated to change the branching ratio by 0.5%. 

The decay was first observed by the CERN NA31 experiment [[| based on 2 observed 
events and has been confirmed by later measurements [T^ . 



Here we report the preliminary result obtained from the 1999 data taking of the CERN 
experiment NA48. 



II. EXPERIMENTAL SETUP AND DATA TAKING 

This measurement was carried out as part of the CERN experiment NA48 at the CERN 



SPS, which has previously reported measurements of the related decays Kl — > e + e~7 JTT 
and K L —> / u + /i _ 7 |L2]. A detailed and comprehensive description of the detector is in 
preparation [T3| . 



The NA48 experiment is designed specifically to measure the direct CP violation para- 
meter 3?(e'/ e) using simultaneous beams of Kl and K$. To produce the Kl beam, 450 GeV/ c 
protons are extracted from the accelerator during 2.4 s every 14.4 s and 1.1 x 10 12 of these are 
delivered to a beryllium target. Using dipole magnets to sweep away charged particles and 
collimators to define a narrow beam, a neutral beam of 2 x 10 7 Kl per burst and divergence 
±0.15 mrad enters the decay region. The fiducial volume begins 126 m downstream from 
the target and is contained in an evacuated cylindrical steel vessel 89 m long and 2.4 m 
in maximum diameter. The vessel is terminated at the downstream end by a Kevlon-fiber 
composite window and followed immediately by the main NA48 detector. The sub-detectors 
which are used in the Kl — ► e + e~e + e - analysis are described below in order of succession 
(see Fig. 0). 
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FIG. 1. Side-view of the NA48 detector. 



A magnetic spectrometer consisting of a dipole magnet is preceded and followed by two 
sets of drift chambers. The drift chambers are each comprised of eight planes of sense wires, 
two horizontal, two vertical and two along each of the 45° directions. Only the vertical and 
horizontal planes are instrumented in the third chamber. The volume between the chambers 
is filled with helium at atmospheric pressure. The momentum resolution is Ap/p = 0.65% 
at 45 GeV/c. 

Two segmented plastic scintillator hodoscope planes are placed after the helium tank 
and provide signals for the trigger. 

A liquid krypton filled calorimeter (LKr) is used for measuring the energy, position and 
timing of electromagnetic showers. Spacial and timing resolutions of better than 1.3 mm and 
300 ps, respectively, have been achieved for energies above 20 GeV. The energy resolution 



is 



a(E) 



0.035 



24^ ® 0.OO6, with E measured in GeV. 



E y/E ^ E 

A hadron calorimeter composed of 48 steel plates, each 24 mm thick, interleaved with 
scintillator is used in trigger formation and particle detection studies. 



III. TRIGGER AND DATA PROCESSING 

Candidate events were selected by a two-stage trigger. At the first level, a trigger re- 
quiring adjacent hits in the hodoscope is put in coincidence with a total energy condition 
(> 35 GeV), defined by adding the energy deposited in the hadronic calorimeter with that 
seen by the trigger in the LKr calorimeter. The second level trigger uses information from 
the drift chambers to reconstruct tracks and invariant masses. For the 4-track part of the 
trigger, the number of clustered hits in each of the first, second, and fourth drift chamber 
had to be between 3 and 7. All possible 2-track vertices were calculated online. At least 
two vertices within 6 m of each other in the axial direction had to be found. 
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For the determination of the 4-track trigger efficiency, downscaled events that passed the 
first level were used. Alternatively, events triggered with the neutral trigger - based on the 
data of the LKr calorimeter - were selected. The neutral trigger applied the following cuts 
to the events online: < 5 peaks in each projection, total energy > 50 GeV, first moment of 
cluster energies ('center of gravity') < 15 cm, and lifetime uncorrected for deflection in the 
magnet < 4.5 ■ tk s - 

During the experimental runs, roughly 100 Terabytes of raw data with typically 20 kbytes 
per event were recorded. Reconstructed output was stored in a compressed data format, 45 
times smaller. In addition, several streams of data were formed for events accepted by about 
40 filter modules for the analysis of neutral and charged two-pion decays, rare decays, and 
events for the detector calibration. 



IV. DATA ANALYSIS 



The data sample which yields Kl — > e + e e + e also has been used to select Kl — > 
^'^Daiitz and k l -> ^Daiitz^Daiitz normalization events, with n° Dalitz -> e+e"7. The 
vertex was reconstructed from the 4 tracks by requiring that the sum of the squared trans- 
verse distances from the transverse vertex position, weighted by the inverse track momentum, 
be minimal. 

Events were preselected by requiring two positive and two negative tracks with distance 
of closest approach to the vertex < 15 cm. 

All clusters in the LKr were required to be in a fiducial area given by an octogon about 5 
cm smaller than the outside perimeter of the calorimeter and an inner radius of 15 cm. The 
distance to dead calorimeter cells had to exceed 2 cm to ensure negligible energy loss. The 
separation between cluster centers was required to be > 3 cm. The energy of each cluster 
was required to exceed 2 GeV, well above the detector noise of 0.11 GeV per cluster. 

Electron candidates were identified by requiring that cluster centers in the LKr be within 
1.5 cm of the calculated shower maxima based upon the extrapolation of each track. The 
efficiency of this procedure was measured to be (99.7 ± 0.1)% To reject pion showers, 
the ratio of cluster energy to track momentum E/p was required to lie between 0.9 and 1.2. 
The efficiency of this selection was determined to be > 95% |IT| . Those track-associated 
clusters with < E/p < 0.8 were classified as pions. From a study of Ks — > n + Ti~ decays, 
the probability of pions to be wrongly classified as electrons is estimated to be 0.9%; the 
pion classification is passed by 97.5% of all pions. 

The fiducial volume was defined by the axial postion 7.50 m < z vertex < 90 m downstream 
of the Ks target. Within this volume, 4-track vertices were determined with a typical 
axial resolution of 0.5 m, as estimated by the MC. The Kl energy had to be in the range 
50 GeV-200 GeV. 



A. Background Rejection and Selection of — > e + e e + e Candidates 

Candidate events for the decay Kl — > e + e~e + e~ with all tracks identified as electrons 
were selected. In principle, the following four classes of background sources are relevant: 
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• Events with two decays Kl — > 7rez/ occuring at the same time and for which the pions 
were misidentified as electrons. Being due to two coincident kaon decays the invariant 
mass of the system can be around and above the nominal Kl mass. These events are 
largely rejected by requiring a good vertex quality: the closest distance of approach 
of each track w.r.t. the reconstructed vertex had to be < 5 cm. Events with separate 
vertices do not pass the level 2 trigger (see above). 

Finally it was required that the measured times for clusters associated to the electrons 
had to be consistent within 3 ns with each other. A study of sidebands in this time 
distribution shows that the background from this source is negligible. 

• Events Kl —> 7r°7r 7r , where the 7r°'s undergo single or double Dalitz decays or pho- 
tons convert in the material of the detector, so that 2 positive and 2 negative electrons 
are detected. Due to the missing photons, the invariant mass of the e + e~e + e~ sys- 
tem is below the nominal Kl mass. The loose requirement that the square of the 
reconstructed transverse momentum p\ of the reconstructed kaon with respect to the 
line joining the decay vertex and the Kl target must be less than 0.0005 (GeV/c) 2 , 
restricts this background to less than 2.2% . The position of the cut is indicated in 
Fig. 0. The Monte Carlo simulation indicates that 1.5% of the signal events are lost 
by the p\ cut. 
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FIG. 2. Correlation of e + e~e + e~ invariant mass with the squared transverse momentum pj, of 
the reconstructed kaon. 



• Events Kl — > 77 and Kl — > e + e~7, with conversion of the photons in the material 
upstream of the spectrometer also yield invariant masses around the nominal mass. 
Each pair of oppositely charged tracks was therefore required to be separated by > 2 cm 
in the first drift chamber as indicated in Fig. @(a). Note that the conversion probability 



6 



in the material of the NA48 detector is of similar magnitude as that for internal photon 
conversion to a e + e~ pair. As the angular opening of oppositely charged tracks peaks 
mostly at small angles, 60% of the signal events are lost by this cut; according to the 
MC, there is no remaining background with converted 7's. 

Events Kl — > 7r + 7r~e + e~ |TJ|], with the pions misidentified as electrons. Due to the 
low misidentification probability of 0.9% this background is found to be negligible. 



(a) Signal 
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(b) Normalization 
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FIG. 3. Distance d of each pair of oppositely charged tracks at the position of the first drift 
chamber for (a) Kl — ► e + e~e + e~ and (b) Kl — > tt + it~ K Dalitz . Dots with error bars show the data, 
the histogram is the Monte Carlo prediction normalized to the number of observed signal events 
with d > 2 cm. 



The invariant e + e~e + e~ mass plot resulting from this selection is shown in Fig. f|. Note 
the slightly asymmetric shape of the Kl mass peak, which is due to photons radiated off 
the electrons in the final state. 

Finally, a mass window of 0.475 GeV/c 2 < m(e + e~e + e~) < 0.520 GeV/c 2 was set to 
define the final sample. In total, 132 candidate events were selected in this way. From the 
number of events observed below the Kl mass peak and their distribution we estimate that 
the background contribution to the signal region is negligible. 
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FIG. 4. Invariant mass of the e + e~e + e~ system. 
B. Normalization 

The four-track decay K L — > 7r + 7r~7r^ aZii2 , with vr^, aiiiz — > e + e~7, was used for norma- 
lization. Only events that passed the pretrigger discussed above, downscaled by 60, were 
selected. This sample was used to determine the efficiency of the 4-track trigger to be 
95.1 ± 0.2%. The 4-track trigger is required for signal and normalization modes. 

Since both the signal and the normalization modes consist of 4-track events, uncertainties 
due to tracking tend to cancel in the ratio of acceptances. Selection criteria similar to those 
used in the signal mode were applied. In addition, at least one extra cluster in the calorimeter 
not associated with a charged track was required. The invariant mass of the e~e~7 system 
had to be in the range of 0.115 - 0.155 GeV/c 2 . Monte Carlo studies showed that 66 % 
of the reconstructed Kl — > 7r + 7r~7r£, candidate decays are from Kl — > 7r + 7r~7r° with one of 
the external photons converting in the material of the detector (see Fig. |](b)). All other 
backgrounds have been estimated to be negligible. Using these cuts, 17123 Kl — > 7r 7r°7r£, 
decays were selected. 

In a second analysis, Kl — > ^^Daiitz^Daiitz even ts were selected, yielding 5167 events for 
normalization. While having a complicated topology of eight clusters, these events have the 
advantage that all decay products interact electromagnetically in the detector and that the 
radiative corrections should be similar to those in the signal mode. 

C. Acceptance Determination and Kaon Flux 

For the simulation of the Kl —> e + e~e + e~ acceptance, the matrix element was taken 
from Ref. || neglecting the interference of the two virtual photons. The distribution of 
the angle spanned by the decay planes of the two e + e~ pairs corresponds to a Kl which is 
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assumed to be entirely CP= —1. The PHOTOS [Tj| package has been used to simulate final 
state radiation both for the signal and normalization channels. 

The acceptance for K L — > e + e~e + e~ is calculated to be 7.8% for events generated in the 
range 45 GeV < E Kl < 215 GeV and 5 m < z vertex < 91 m. The normalization to the K L 
signal has been measured from the number of — » ^ + ^~^D a ut z decays in the same sample. 
Using the acceptance of 1.34%, calculated by Monte Carlo simulation, a total number of 
5.1 x 10 10 Kl decays is obtained. With the cuts described above, the inclusion of radiative 
corrections decreased the acceptance of signal and normalization channel by 8.8% and 3.4%, 
respectively. 

Consistent results for the total number of were obtained when the alternative nor- 
malization channel, Kl — > ^^Daiitz^Daiitzy was used instead. 



V. RESULTS AND DISCUSSION 

A study of the stability of the branching ratio determination was made as a function 
of the cuts applied. A systematic error of 3.5% was estimated, mainly being due to the 
variation in the minimal distance of clusters from the beam pipe and the cut on the vertex 
quality. 

A second contribution comes from the 4-track trigger inefficiency. A Monte Carlo si- 
mulation of the level 2 algorithm yields a 95.0% efficiency for the normalization mode, in 
good agreement with the measured value. For the signal, the simulated efficiency is higher 
(99.8%). We chose not to apply a correction to the branching ratio; instead we introduced 
a systematic error of ±5% for this preliminary result. 

Finally, the effect of overflows in the drift chambers has been considered. If events with 
an overflow condition in a window of 312 ns around the event time are removed, 20% of the 
events are lost and the branching ratio stays constant within 1%. 

Adding these sources of systematic error in quadrature, we obtain a total systematic 
error of ± 6.2%. 

From the numbers given above, a branching ratio of 

is obtained, where the statistical and systematical uncertainties as well as the uncertainty 
in the branching ratio of the normalization channel are given separately. 

This result is consistent with the theoretical expectation of || and the previous aver- 
age value of (4.1 ± 0.8) x 10~ 8 flTSI , with 5 times the statistics of the single best previous 
experiment. 
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